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Abstract 

Data on macroscopic superconducting anisotropy of MgE>2 are reviewed. The data 
are described within a weak coupling two-gaps anisotropic s-wave model of super- 
conductivity. The calculated ratio of the upper critical fields 7^ = H C 2 A bl 'Hc2,c 
increases with decreasing temperature in agreement with available data, whereas 
the calculated ratio of London penetration depths 7^ = A c /A a ;, decreases to reach 
1.1 at T = 0. Possible macroscopic consequences of 7a 7^ 7# are discussed. 
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1 Introduction 



Physics and applications of new superconducting materials cannot be prop- 
erly understood and developed without careful characterization of possible 
anisotropics. High-T c materials are one of the best examples of relevance of 
the superconducting anisotropy to all aspects of physics and possible applica- 
tions. 

It is a common practice to characterize anisotropic superconductors by a single 
anisotropy parameter denned as 7 = £ a /£ c = A c /A a . Here, £ is the coherence 
length, A is the penetration depth, and a, c are principal directions of a uni- 
axial crystal of the interest here. Often £ a = are denoted in literature as 
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£ab] for formal reasons we prefer a single subscript notation. Since the up- 
per critical fields, H c2c = 0o/2vr^ and H c2a = 0o/27r£ a £c, the anisotropy 
parameter can also be written as 7 = H c2>a / H c2c . Historically, the practice 
to assign a single anisotropy parameter to each material emerged after the 
anisotropic Ginzburg-Landau (GL) equations were proposed phenomenologi- 
cally by Ginzburg [1] and derived microscopically for arbitrary gap and Fermi 
surface anisotropies by Gor'kov and Melik-Barkhudarov [2]. Formally, this 
came out because near the critical temperature T c , the same "mass tensor" 
rriik determines the anisotropy of both £ (of the upper critical fields H c2 ) and 
of A: 

(£ 2 ) ifc oc m~f} , (\ 2 ) ik oc m ik . (1) 



Therefore, 



At arbitrary temperatures, however, the theoretical approach for calculating 
H c2 (the position of the second order phase transition in high fields) has little 
in common with evaluation of A (the weak-field relation between the current 
and the vector potential), so that the anisotropies of these quantities are not 
necessarily the same. In fact, in materials with anisotropic Fermi surfaces 
and anisotropic gaps, not only the parameter 7^ = H c2 a /H c2 c may strongly 
depend on T, but this ratio might differ considerably from 7a = A c / A a at low 
temperatures [3,4]. 

In this brief review we collect data on anisotropy of MgB2 available to us, 
discuss briefly the methods used and the results obtained. Then we outline 
the weak-coupling microscopic approach which can be used to evaluate the 
anisotropy parameters and their T dependence. We conclude with discussion of 
macroscopic consequences of different H c2 and A anisotropies focussing on the 
torque problem, the quantity commonly measured to extract the anisotropy 
parameters. 



2 Data review 



Experimentally, superconducting anisotropy of MgB 2 was under scrutiny right 
away after the discovery of this material. Almost all experimental studies were 
concerned with the H c2 anisotropy, 7^ = H c2 ^ a b/ H c2tC . 

Based on the AC susceptibility and the field dependent magnetization mea- 
surements on sorted powders, magnesium diboride was claimed to be an 
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isotropic superconductor [5]. The H C 2 anisotropy for separate particles set- 
tled on a flat surface was reported as 7# ~ 1.73 [6]. For a hot-pressed bulk 
samples values of up to 1.2 were obtained from susceptibility and resis- 
tance measurements [7]. Magnetotransport measurements on MgB 2 thin films 
with different degrees of c-axis orientation yielded t# — 1-8 — 2 (T c — 31-37 
K) with higher anisotropy for films of higher resistivity [8]. The temperature- 
independent 7# = 1.2 (T c = 39 K) was reported for c-oriented films, Ref. [9]. 
Similar results {^h = 12- 1.8) were later obtained for a set of three films 
pulsed laser deposited (PLD) on different substrates (T c = 31-37 K) [10,11]. 
In this case no apparent correlation was observed between the H c2 anisotropy 
and residual resistance ratio (RRR) or superconducting transition tempera- 
ture of the films. Quite different values, 7h = 9 - 13, were reported for in-situ 
grown PLD films with T c suppressed down to 24 - 27 K [12]. 

In a different approach, the H c2 anisotropy of magnesium diboride was evalu- 
ated from a number of measurements on randomly oriented powders. Analysis 
of the conduction electron spin resonance (CESR) data on fine MgB2 powders 
taken in a wide temperature range yielded an estimate of 7^ ~ 8 [13]. Similar 
idea (deconvoluting the signal in a certain temperature range in two compo- 
nents, one corresponding to superconducting state and the other originating 
from the normal state) was later used for interpretation of the 11 B nuclear 
magnetic resonance (NMR) data in the mixed state of MgB2 [14] . These mea- 
surements resulted in 7^ pa 6. Subsequent n B NMR and magnetization studies 
on Mg 1 _ a; Al x B2 (x < 0.025) reveal a considerable decrease of the anisotropy 
with Al doping [15]. 

In addition to local NMR and CESR probes, global measurements on samples 
with randomly oriented grains were used for the anisotropy studies. High field 
reversible DC magnetization data were well described assuming that MgB 2 is 
an anisotropic superconductor with 7^ = 6 - 9 and using the available ex- 
perimental estimates of the average penetration depth A(0) = (A^ b A c ) 1//3 = 
110 — 140 nm [13]. The H C 2 anisotropy was also evaluated from the broaden- 
ing of the resistive transition in applied magnetic field. Two slightly different 
theoretical approaches to this problem, [16] and [17], based on the anisotropic 
Ginzburg-Landau and percolation theories for randomly oriented anisotropic 
superconducting grains, were developed more than a decade ago in early days 
of the high T c superconductivity. They require somewhat different analysis of 
the resistivity data but result in similar estimates: 7# = 5 - 7 [18]. 

The aforementioned determinations of the H c2 anisotropy utilized traditional 
techniques, with the exception of those based on CESR and NMR. About a 
year ago a robust and a simple way to extract this anisotropy from the data 
on T dependence of the magnetization of random powders taken in various 
fixed fields has been suggested [19]. The method is based on two features in 
(dM / dT) H : the onset of diamagnetism at T™ ax , which is commonly associated 
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with H C 2 and which is, in fact, H C 2 jma x = Hc2,a, and a kink in OM/dT at a 
lower temperature T c mm (see Fig. 1). The origin of these two features can be 
understood as follows. 

Upon cooling a powder sample in a fixed H, there is a deviation from a roughly 
T independent normal state magnetization (for non-magnetic materials as is 
the case of MgE^) to an increasingly diamagnetic signal at T = T™ ax . On the 
H — T diagram, the point (H,T™ ax ) lies at the phase boundary H C 2 tTnax (T). 
Repeating this measurement at various fixed field one recovers the full curve 
H C 2, m ax(T) = H c2: a(T). This line coincides with the upper critical field deter- 
mined for polycrystalline samples by, e.g., the standard resistivity measure- 
ments. A second sharp feature occurs in OM/dT when the sample tempera- 
ture passes through T c mm . This can be readily understood by considering what 
happens to the sample upon warming. For T < T c mm all randomly oriented 
grains are superconducting, whereas for T > T c mm part of them are normal, 
depending upon their orientation with respect to the applied magnetic field. 
Therefore, upon warming through T c mm we expect some peculiarity, or a kink 
in dMj dT associated with the onset of normal state properties in an increasing 
number of appropriately oriented grains. In a similar spirit the field dependent 
magnetization can be analyzed. This method constitutes a robust procedure 
independent either of a particular model for describing the anisotropy or of 
a degree of powder randomness. A kink in dM/dT at T = T c mm should be 
present for any angular distribution of the grains (as long as it is continuous, 
although not necessarily random). 

It should be noted that in superconductors with strong fluctuations (like 
HTSC) or in the samples with a distribution of T c 's or other physical properties 
due to chemical or/and structural inhomogeneities the region of superconductor- 
normal transition and the kink in dM/dT may be smeared. In such cases this 
simple method may not work. Still, given difficulties in growing single crystals 
of a sufficient size often encountered when a new material is synthesized, the 
method may prove useful and sometimes the only one available for determi- 
nation of the H c2 anisotropy. 

The analysis of detailed temperature and field dependent magnetization mea- 
surements on MgB 2 powders together with the data on polycrystalline H c2 
(= H" 2 mx ) obtained from high field magnetotransport measurements on simi- 
lar samples yielded a complete anisotropic H C 2(T) phase diagram [20], see Fig. 
2. 

A number of direct measurements of the H c2 anisotropy on single crystals was 
reported. Magnetoresistance data on small MgB2 single crystalline platelets 
[21,22,23,24] gave ^ H = 2.6 - 3.0 (in these samples T c = 38.0-38.6 K, RRR = 
5 - 8). It is argued in Ref. [25] that electrical transport measurements are not 
well suited to probe bulk upper critical field of MgB 2 due to possible presence 
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of a second phase (related to surface effects) with enhanced superconducting 
parameters. In this paper, the bulk if" 2 (T) and H^ 2 (T) are evaluated utilizing 
data on the in-plane thermal conductivity for two orientations of the applied 
magnetic field. The extrapolated value of 7^ is 7#(T — > 0) = 4.2 [25]. 

The torque magnetometry was also employed to study superconducting anisotropy 
on single crystals [26]. From angular dependent torque measurements taken 
in different applied fields the anisotropy of the upper critical field was found 
to be temperature dependent: it decreases from 7^ ps* 6 at 15 K to ps 2.8 at 35 
K. Reversible angular dependent torque measurements (with a vortex shaking 
process [27]) performed in the temperature range 27 - 36 K and in fields up 
to 10 kOe were analyzed using a formula derived in [28] for 7^ = 7a, see Eq. 
(15) below. As a result, nearly linear field dependence of 7 has been extracted 
with 7(ff — > 0) pa 2 and 7(10 kOe) ps 3.7 with practically no T dependence 
between 27 K and 36 K. Extensive torque data of Ref. [29] were also inter- 
preted assuming 7# = 7a = 7 and produced estimates 7 ~ 3 — 4 which vary 
somewhat with T and H . 

Subsequent analysis [20] of the anisotropic H c2 (T) data obtained on random 
MgB 2 powders (Fig. 2) showed 7#(T) dependence consistent (Fig. 3) with the 
one observed in the single crystals [26] . Similar temperature-dependent be- 
havior was extracted from magnetotransport, magnetization, ac susceptibility 
and specific heat in applied magnetic field measurements on single-crystalline 
platelets [30,31,32] and in thin films [11]. 

Summarizing, the reported anisotropy of the upper critical field in magnesium 
diboride range from 7^ ps 1 to 7^ ps 13. There are some concerns related to 
the spread of these results such as purity and homogeneity of materials used, 
degree of alignment of crystallites/grains, possible sharp angular dependence 
of H c2 so that correct evaluation of the anisotropy may suffer from even slight 
misorientations in measurements on single crystallites, and, last but not least, 
adequacy of models used to analyze the data on polycrystals and crystals. Even 
for similar small single crystals different apparently direct measurements give 
somewhat different estimates of 7^ (see e.g. [26,25]). On the other hand, there 
is an apparent consistence between recent data on the temperature dependence 
of 7h obtained by at least two groups on single crystals and on high quality 
powders, similar "convergence" emerges in data for other physical properties. 
This indicates that we are coming closer to measuring the intrinsic properties 
of MgB 2 . As a side remark it should be mentioned that temperature dependent 
H c2 anisotropy is not unique for MgB 2 , it was observed, for example, in NbSe 2 
[33]. 

For sub- mm size MgB 2 single crystal, an anisotropic lower critical field, H c i, 
was estimated from the M(H) data [21]. The H^(T) dependence was reported 
as near linear with H^(0) ps 384 G, while the H c cl (T) appeared to be non-linear 
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with and H^(0) ~ 272 G. Another estimate of the temperature dependent 
anisotropic H c i from magnetization measurements on MgB 2 single crystallites 
[32] gave close to linear if cl (T) for both directions with extrapolated to T = 
values H^(0) « 2800 G and if#(0) « 1300 G and near T independent ratio 
H°JH* « 2.2. 



Apparently more work on sample quality, measurement procedures, measure- 
ments analysis, and modelling is required to obtain a better physical picture 
of the superconducting anisotropy of MgB2. Reliable data on the tempera- 
ture dependence of the anisotropy of the London penetration depth, when 
available, will serve as an important test of existing models. 



3 Theory of anisotropies of A and H c2 



To address theoretically the question of anisotropy at arbitrary temperatures 
one has to turn to a microscopic model. The most complete and sophisticated 
approach is based on taking into account the actual Fermi surface, the phonon 
distribution, and the electron-phonon interaction (all of which are anisotropic) 
in the frame of the Eliashberg theory. This can be and has been done numeri- 
cally by a few groups, which evaluated the Fermi surfaces [34,35] and the gap 
distribution [36,37]. 

As is shown in Fig. 4, the Fermi surface of MgB 2 consists of two nearly cylindri- 
cal sheets, "2D a-bands", and two others corresponding to the "3D 7r-bands". 
The gap on the four Fermi surface sheets of this material has two sharp max- 
ima: Ai ~ 1.7 meV at the two 7T-bands and A 2 ~ 7meV at the two cr-bands, 
see Fig. 5. Within each of these groups, the spread of the gap values is small, 
and the gaps can be considered as constants, the ratio of which is nearly T 
independent. In this situation, a model with only two Fermi surface sheets and 
two gaps may prove useful in relating various macroscopic properties of MgB 2 . 
Starting with Ref. [38], the two-band models were studied by many, see, e.g., 
Ref. [39] and references therein. We outline below the results of evaluation of 
the London penetration depth and the upper critical fields within a relatively 
simple two-gap weak-coupling model which provides consistent physical pic- 
ture of superconducting anisotropies. We also discuss possible consequences 
of different anisotropies of A and H c2 upon macroscopic magnetic properties 
of MgB 2 . 
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3.1 Two gap model 



We begin with the notion that within the weak-coupling BSC theory the 
macroscopic characteristics of superconductors such as the coherence length 
and the penetration depth are expressed in terms of quantities averaged over 
the Fermi surface. This basic consequence of the collective nature of the super- 
conducting condensate is formally expressed by the self-consistency (or gap) 
equation: 



Here, N(0) is the total density of states at the Fermi level per one spin; the 
Fermi velocity v relates to a certain position on the Fermi surface, A is the 
gap function, f(r,v,u) is Eilenberger Green's function which describes the 
condensate, and hoo = nT(2n+ 1) with an integer n. In the uniform zero-field 
situation / = A/VA 2 + h 2 uo 2 . Substituting this in Eq. (3) one can evaluate 
how the gap depends on temperature and on position at the Fermi surface 
provided the pair-pair interaction V(v, v') is known. 

The averages over the Fermi surface weighted with the local density of states 
oc 1/v are defined as 



Commonly, the interaction V is assumed factorizable [40], V(v, v') = V f2(v) Q( 
(the assumption with far reaching consequences for the two-gap model, see be- 
low). Then, one looks for A(r, T; v) = \P(r, T) Q(v). The self-consistency Eq. 
(3) for the uniform condensate takes the form: 



A(r, v, T) = 2nTN(0) £ (V(v, v')f(r, v' , u))* . 



(3) 




(4) 



1 = 2irTN(0)V £ ( 




(5) 



Consider now a model material with the gap anisotropy given by 



fi(v) = fti • 



v6f] 



1,2 ) 



(6) 



where F 1 ,F 2 are two sheets of the Fermi surface. Denoting the densities of 
states on the two parts as Aq i2 , and assuming the quantity X being constant 



at each sheet, we obtain for the general averaging: 



(XjiNi + (X) 2 N 2 
N(0) 



v 1 (X) 1 + v 2 (X) 



2 ■ 



(7) 



where we use normalized densities of state u^ 2 = Ni i2 /N(0) for brevity. The 
Fermi relief of the gap, fi, can be normalized so that (f2 2 ) = 1 which gives for 
the two-gap model: 



Within this model, the interaction V now takes the form of a 2 x 2 symmetric 
matrix Vij/Vo with diagonal components Qf and Q 2 , and the two equal off- 
diagonal elements QiQ 2 . It is worth noting that if V i2 = 0, i.e., there is no 
interband interaction at all, the general gap equation (3) splits in two indepen- 
dent equations for the independent order parameters, A x on Fi and A 2 on F 2 . 
Such a model would imply different, in general, critical temperatures for A's, 
different phases with many macroscopic consequences [41,42]. The factoraz- 
able interaction Vij = VoQiQj excludes such possibilities and automatically 
implies the single T c for both gaps, the feature in fact seen in many tunneling 
experiments on MgB 2 [43,44,45]. 

3.2 T dependence of the two gaps 

The self consistency equation (5) for the order parameter ^(T) can be solved 
for known f2i j2 and the densities of states v\^ 2 . Based on the band structure 
calculations, the relative densities of states v\ and u 2 of our model are ~0.56 
and 0.44 [46,36]. The ratio A 2 /Ai = f2 2 /f2i ps 4, if one takes the averages of 
6.8 and 1.7 meV for the two groups of distributed gaps as calculated in Ref. 
[36]. Then, the normalization (8) yields Q,\ = 0.36 and Q 2 = 1.45. It should be 
noted that the tunneling data provide somewhat lower ratio A 2 /Ai [43,44,45]. 

We have now all parameters needed to solve the self-consistency equation (5) 
for \&(T) (the clean case). This is done numerically and the result is shown in 
Fig. 6 along with two gaps A lj2 (T). 

The parameter A/T c commonly used to distinguish strong and weak coupling 
materials, has a "strong coupling" value of 4 for A 2 /T c , whereas Ai/T c m 1, 
which is less than the BCS weak coupling value of 3.5 . Thus, this parameter 
cannot be used to characterize coupling for materials like MgB 2 . The T de- 
pendence of the two gaps similar to shown in Fig. 8 has been reported by a 
number of groups [43,44,45]. 



Vl\vi + Vl\v 2 = 1 , vi + v 2 — 1 . 



(8) 



8 



3.3 Anisotropy of A 



Having obtained the temperature dependence of both gaps, one can turn to 
calculation of the weak- field penetration depth A(T). To this end, one first 
solves equations of superconductivity for the zero-field case and then employs 
the perturbation theory to evaluate the small corrections due to small currents. 
In this manner one obtains the London relation between the current and the 
vector potential, 

4 71 " . i /27T 89 . \ , n . 



where 9 is the phase, 0o is the flux quantum, A is the vector potential, and 



is the tensor of the inverse squared penetration depth. Here, i,k — x,y, z and 
f3 2 = A 2 + h 2 uj 2 (the reader may find details of this calculation in Ref. [3]). 
This general expression for the tensor of squared penetration depth is valid for 
clean superconductors with an arbitrary Fermi surface and gap anisotropics. 

The anisotropy parameter now reads: 

2 _ Ag c _ (^AgE,r 3 ) (u) 
7A Ka (^A 2 E./3- 3 )' 1 ) 



As T -> 0, we have 2ttTA 2 £ w (3~ 3 -> 1, and 

7^(0) = kf. (12) 



c / 



Note that the gap and its anisotropy do not enter this result. The physical 
reason for this is in the Galilean invariance of the superfluid flow in the absence 
of scattering: all charged particles take part in the supercurrent. 

Near T c , the sum over uj in Eq. (11) is is a kp independent constant because 
P ftw, and we obtain the GL result of Ref. [2]: 

^) = AL = (?S|)' < 13 > 

We observe that near T c , the gap anisotropy amplifies contribution of the large 
gap Fermi surface pieces to the anisotropy parameter 7 A . We also see that for 
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isotropic gaps, jx(T c ) = 7a (0) (which is not true for the general case). Thus, 
the anisotropy parameter depends on T, the feature absent in superconductors 
with isotropic gaps. 

To apply these results to MgB 2 , one has to know averages of squared Fermi 
velocities not only over the whole Fermi surface (as for 7a (0) of Eq. (12)), but 
also the averages over the separate Fermi sheets. These can be taken from the 
band structure calculations of Ref. [46]: 

(Oi = 33.2, (v 2 ) 1 = 42.2, 

{v 2 a ) 2 = 23 , (v 2 c ) 2 = 0.5 x 10 14 cm 2 js 2 . (14) 
The numerical results for 7a (T) are shown in Fig. 7. 

Thus, we expect the London penetration depth of clean MgB 2 to be nearly 
isotropic at low temperatures; J\(T) increases on heating and reaches ~ 2.6 at 
T c . Qualitatively similar behavior of J\(T) is predicted within a more general 
Eliashberg scheme [47]. As up to date, there is no published data on direct 
measurements of the A anisotropy and its T dependence. It is pointed out 
below that the data on the angular dependence of the torque acting on the 
anisotropic superconducting crystal in tilted fields cannot provide information 
about jx(T) without taking into account the H c2 anisotropy which for MgB 2 
might be quite different from jx(T). In principle, one can extract 7a from the 
data on the vortex lattice structure in fields tilted relative to the c axis of 
MgB 2 . These, however, are difficult to acquire and to our knowledge are not 
yet available. 



3.4 Anisotropy o/H c2 

Unlike the case of London A, the evaluation of H c2 is more complicated a 
problem because one has to deal with the high-field phase transition from the 
normal to the superconducting mixed state. Still, major behavior of H c2 (T) 
can be obtained within a simplified scheme: Since the actual band structure 
enters equations for H c2 via Fermi-surface averages, one can further model the 
sheets F 1>2 by two spheroids with average Fermi velocities close to the band- 
structure generated values. As a result one obtains qualitatively (and - given 
the spread of existing data - quantitatively) correct behaviors of H c2 (T) for 
both principal directions. Details of this evaluation can be found in Ref. [4]. 
The result for the anisotropy parameter 7#(T) = H c2>a b/H c2)C obtained within 
the scenario of two constant gaps on two Fermi surface sheets, the same used 
for 7a (T) and with the same input parameters, is shown in Fig. 7. 
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The drop of the H C 2 anisotropy with increasing T has been recorded in mea- 
surements on single crystals of MgB 2 [26] . Nearly the same behavior has been 
extracted from the data on random polycrystals in the whole temperature 
range using the T dependence of the magnetization of random powders [20]. 
Recent specific heat measurements [31] and magnetization data [30] on sin- 
gle crystals produced similar results. All these data show qualitatively similar 
behavior to that of the upper curve of Fig. 7. 

Physically, the large anisotropy of H c2 at low temperatures (~ 6 in our cal- 
culation) is related to the large gap value at the Fermi sheet which is nearly 
two-dimensional. With increasing T, the thermal mixing with the small-gap 
states on the three-dimensional Fermi sheet suppresses the anisotropy down 
to 2.6 at T c . 

3. 5 Free energy and torque in materials with different anisotropies of A and 



One of the most sensitive methods used to extract the anisotropy parameter 
7 is to measure the torque acting on a superconducting crystal in the mixed 
state with the applied field tilted relative to the crystal axes. In intermediate 
field domain, H c i <C H <C H c2 , the demagnetization shape effects are weak, 
and the torque density can be evaluated: [28,48] 



where 6 is the angle between the induction B and the crystal axis c, 



A 3 = A^A C , and r\ ~ 1. This formula can be written as r = M x H; since in 
this field domain the magnetization M <C H , one can disregard the distinction 
between B and H. It has been assumed in the derivation of Eq. (15) that 
the anisotropies of and of the London penetration depth are the same: 
1h = 7a = 7- It is worth noting that Eq. (15) describes the system with the 
stable equilibrium at 9 — ir/2, i.e., the uniaxial crystal in the external field 
positions itself so that the field is parallel to ah (if 7 > 1). 

Expression (15) for the torque has been derived within the London approach by 
employing the cutoff at distances on the order of the coherence length £ where 
this approach fails. The formula, however, has been confirmed experimentally 
(for a few high-T c materials close to their T c ) with a good accuracy as far as the 



0off( 7 2 -l)sin2fl rjH^a 
647r 2 A 2 7 1 /3 £ (#) 11 Bg (0) > 



(15) 




(16) 
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angular dependencies of quantities involved are concerned [49]. Uncertainties 
of the London approach are incorporated in the parameter r\ ~ 1; discussion 
of those can be found, e.g., in Ref. [50]. 

We now outline possible effects of different anisotropies of A and H c2 on the 
torque angular dependence. To this end, we first recall the London expression 
for the free energy valid for intermediate fields H c \ <C H <C along z tilted 
with respect to the c crystal axis over the angle 9 toward the crystal direction 
a [51]: 

F _ B 2 | &m zz y 1 {m 

8tt 8vrA 2 m a ^ m zz G 2 x + m c G 2 y ' 1 ' 



here, m zz = m a sin 2 9+m c cos 2 9, m c /m a = 7 2 , m 2 a m c = 1 for uniaxial crystals, 
and G form the reciprocal vortex lattice. The summation is extended over all 
nonzero G. 

For isotropic case, all reduced masses in Eq.(17) are unity. The sum G~ 2 
can be approximated by an integral (<po/47r 2 B) J dG x dG y /G 2 which diverges 
logarithmically and should be cut off at a circle in the reciprocal space of a 
radius ~ l/£ where £ is the radius of the circular vortex core. 

In anisotropic material with 7a = 7^, calculation of this sum is reduced to 
the isotropic case by rescaling 

gx = y/mZG x , g y = y/rr^G y . (18) 



In this manner one can derive the free energy of a superconductor in an applied 
field of a certain direction; the torque formula (15) is obtain differentiating 
the energy with respect to the angle 9. 

However, for 7a 7^ 7h, the masses in the sum (17) correspond to the anisotropy 
7a (in other words, the vortex-vortex interaction is determined by the A 
anisotropy) whereas the cut-off reflects the shape of the core (it is determined 
by the anisotropy of the coherence length £, i.e., by the parameter 7^). In this 
situation, the scaling (18) which transforms the summand of Eq. (17) to the 
isotropic-like form, leaves the shape of the cut-off contour elliptical. Therefore, 
the energy F and its angular dependence acquire extra terms and so does the 
torque which can be measured. 

We refer the reader to the original work [52] for the details of this calculation 
and give here only the result. To write explicitly the angular dependence of 
F, it is convenient to use the angular functions 

®x,h(9) = ex, H {9)hx, H (19) 
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where 6x^(9) are denned in Eq. (16) with corresponding 7's. In terms of these 
functions, 

_ B 2 <t> BQ x 2v^/x a 0oe A 

" 8vr + 32^ (6 A + Qn? ' ( j 



where \i a = 7# 2 ^ 3 . The torque density r = —dF/d6 follows: 
foB^l-l) sin 2fl r ^ff c2 , c 4e 2 6 A x 



T = 



where e = 2.718... . In the standard case of = lx = 7, Eq- (21) reduces to 
the result (15). 

Since both A # < 0, the second contribution to the torque (21) is negative 
whereas the first one is positive. The positive torque implies that the system 
energy decreases with increasing 9, as in the standard case of 7^ = 7 A for 
which 9 = 7r/2 is the stable equilibrium. 

On the other hand, it is straightforward to see that for large enough 7^ and 
fixed 7 A (e.g., set 7 A = 1 and leave 7# > 1) Eq. (21) yields a negative torque, 
which means that in the equilibrium the crystal c axis is parallel to the applied 
field. This strange behavior can be tracked down to the angular dependence 
of the line energy of a single vortex: the energy cost of having an elliptical 
core for vortices along ab for which the current distribution far from the core 
is near isotropic (as for the example with 7 A = 1) is too high. 

To illustrate how the angular dependence of the torque varies with anisotropies 
of H c2 and A, we evaluate numerically the torque density (21) for parameters in 
the range of those for clean MgB 2 . Fig. 8 shows t{9) for 7 A = 2.2 and 7# = 3, 
the values expected for temperatures somewhat below T c . Qualitatively, the 
dependence is standard; the torque is positive in the whole angular domain, 
i.e., 6 = 7r/2 is the stable equilibrium. 

With decreasing T, 7^ of MgB 2 increases whereas 7 A decreases. In Fig. 9 the 
torque (21) is plotted for 7 A = 2, ^ H = 5 (the upper curve) and for 7 A = 1.7, 
7// = 5.3 (the lower curve). These values correspond roughly to 0.7 and 0.6 T c 
according to Ref. [4]. Clearly, 9 = n/2 as well as 9 = are unstable; the stable 
equilibrium is shifted to0<#<7r/2. 

Finally, we plot in Fig. 10 the torque density for parameters which correspond 
to low temperatures, where 7^ « 1.1 and 7^ ~ 6. The torque is negative for 
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all angles implying the stable equilibrium at 9 = 0. 



4 Conclusions 



Hence, different gaps at different Fermi surface pieces of MgB 2 (or gener- 
ally, anisotropic gaps on anisotropic Fermi surfaces) may lead to profound 
macroscopic consequences such as those considered above. Moreover, since 7# 
determines the anisotropy of the coherence length and, therefore, of the vortex 
core, whereas 7a describes the ellipticity of the current distribution far from 
the core, the symmetry of the intervortex interaction should depend on the 
intervortex spacing, i.e., on the field B and its direction. In clean MgB 2 at 
low temperatures in fields along ab, one expects to have the standard trian- 
gular (hexagonal) vortex lattice in low fields (7a ~ 1), which should evolve 
to a distorted triangular (orthorhombic) lattice in increasing fields when the 
core ellipticity (7^ ~ 6) becomes relevant for the vortex current distribution. 
The field dependence of the vortex lattice structure should become weaker at 
elevated temperatures and disappear near T c where 7a — > 7#. In intermedi- 
ate temperature range where the equilibrium vortex lattice orientation shifts 
from 9 = 7r/2 to lower angles, one may expect peculiar vortex dynamics (as 
for fields near parallel to the layers of high-T c materials). This possibility calls 
for further study. 

Thus, the question "what is the anisotropy parameter of MgB 2 ?" does not 
have a unique answer. To pose the question properly one should specify the 
quantity of interest. If this is H c2 , the answer is given by the upper curve 
of Fig. 7 for a clean material; if this is A, see the dashed line. If this is the 
magnetization in intermediate fields, M oc (0o/A 2 ) ln(if c2 /if), the main con- 
tribution to anisotropy comes from A; however, the if c2 anisotropy contributes 
as well (being smoothed by the logarithm). The last situation should be taken 
into account while extracting anisotropy from the torque data in tilted fields 
[26,29]. It should be noted in closing, that all anisotropies discussed here might 
be suppressed by impurities. 
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Fig. 1. Sketch of the maximum H C 2 )m ax{T) = H C 2,a an d the minimum 
H C 2,min(T) = H C 2 )C upper critical fields. For a given applied field H, the relation 
H = # C 2,mm(r c mi ") = H c2:max (T™ ax ) defines temperatures T™ in ,T™ ax . The open 
(shaded) rectangles represent superconducting (normal) grains. 
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Fig. 2. The H — T phase diagram for anisotropic MgE>2. Symbols: circles (open 
and filled) are extracted from M(T)\h, triangles - from M(H)\t, asterisks - from 
resistivity data on polycrystalline MgB2- 
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Fig. 3. Temperature-dependent anisotropy of the upper critical field. The range of 
data from Ref. [26] is shown for comparison as a hatched area between dashed lines. 
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Fig. 5. The energy distribution of the gap values calculated in Ref. [36]. 
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Fig. 10. The same as in Fig.8, but 7^ = 1.1 and 7^ = 6. 
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